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Abstract 

We introduce a new definition of 7r-flatness for linear differential delay systems 
with time-varying coefficients. We characterize ir- and 7r-0-flat outputs and provide 
an algorithm to efficiently compute such outputs. We present an academic example of 
motion planning to discuss the pertinence of the approach. 
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1 Introduction 

Differential flatness, roughly speaking, means that all the variables of an under-determined 
system of ordinary differential equations can be expressed as functions of a particular output, 
called flat output, and a finite number of its successive time derivatives ([BUI [3], see also 
[U O |S] and the references therein) . 

For time-delay systems and more general classes of infinite-dimensional systems, exten- 
sions of this concept have been proposed and thoroughly discussed in [TJ [HI El EI] ■ In a linear 
context, relations with the notion of system parameterization [11[I12| and, in the behavioral 
approach of [13], with latent variables of observable image representations [M], have been 
established. Other theoretic approaches have been proposed e.g. in [TSJ HE]- Interesting 
control applications of linear differential-delay systems may be found in pj [TU] . 
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Characterizing differential flatness and flat outputs has been an active topic since the 
beginning of this theory. The interested reader may find a historical perspective of this ques- 
tion in 0|6]. Constructive algorithms, relying on standard computer algebra environments, 
may be found e.g. in [TTJUH] for nonlinear finite-dimensional systems, or [19, 20J for linear 
systems over Ore algebras. 

For algebraic reasons recalled in Section [4j the notions of 7r-freeness and 7r-flatness have 
been introduced in the context of linear differential-delay systems by El EH and 
then for linear time-varying differential-delay systems, with polynomial time dependence, 

byrjBiEg. 

The results and algorithms proposed in this paper concern time- varying differential-delay 
systems with meromorphic time-dependence. Our approach may be seen as an extension of 
[2H EH 03 to this context. 

Our main contributions are (1) a new definition of n- and 7r-/c-flatness, (2) the charac- 
terization of 7r- and 7r-0-flatness in terms of the hyper-regularity of the system matrices, (3) 
yielding an elementary algorithm, of polynomial time complexity, to compute n- and 7r-0-flat 
outputs, based on the row/column reduction of the former matrices. 

Let us emphasize that the introduced definitions are directly suitable for motion planning 
and that our formalism is able to consider state and input variables separately or not. Note 
that, though the natural algebraic framework used to describe such systems is based on a 
ring, denoted by O, of polynomials in both the delay and differential operators, which is non 
principal, the evaluation of our n- and 7r-0-flatness criteria relies upon computations over 
a principal ideal ring, denoted by 0, containing 0, thus making the computations much 
simpler. 

The paper is organized as follows. A first section briefly presents recalls on operators 
and signals (Section [2]), followed by recalls on matrices over the non commutative differential 
ring O (Section[3]). In the latter section, the notions of hyper-regularity and row and column 
reduction are introduced and characterized in this non commutative context. Section [4] then 
deals with differential algebraic notions of systems, 7r-flatness and 7r-fc-flatness and contains 
the main results characterizing 7r-flat and 7r-0-flat outputs. Algorithms to compute such ir- 
flat and 7r-0-flat outputs are then deduced. Finally, the proposed methodology is illustrated 
by an example of motion planning in Section [5j 



2 Operators Sz signal space 
2.1 Recalls on Ore polynomials 

In this paper we consider linear mixed differential and time-delay operators. We will model 
them using the so-called Ore polynomials. These are a class of non-commutative polynomi- 
als, named after Oystein Ore who was the first to discuss them in |24j . In the following we 
will give a brief introduction to Ore polynomials, pointing to the appropriate references for 
details. 

Let K be a ring and let a : K — > K be an automorphism. An additive map i? : K — > K is 
called a a -derivation if for all a and b € K the a-Leibniz rule ??(a&) = cr(a)$(6) +{)(a)b holds 
(compare with [25] Sect. 7.3]). Consider the free If- left module generated by the powers of 
an indeterminate x. We define the right-multiplication of x by an element of K with the 
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commutation rule 

xa = a(a)x + i?(a) for all a £ K. 

Assuming associativity and distributivity, this rule allows us to compute arbitrary products. 
It can be shown (see, e.g., [551 Thm. 7.3.1]) that this makes the free module into a ring 
which we call the ring of (left) Ore polynomials in x w.r.t. a and In the literature this 
is usually denoted by K[x; a, i?]. 

The degree of an Ore polynomial p is defined as the largest exponent n such that x n has 
a non-zero coefficient in p. We use deg = — oo. If K is a domain, then we have the familiar 
rule deg(pq) = degp + deg q for all p and q £ K[x; a, If if is a division ring, then it is 
possible to divide elements in K[x; a, i?] with remainder in a way which is very similar to the 
usual polynomial division. See [53] for details. This turns K[x\a, #] into a (left and right) 
principal ideal domain and thus into a (left or right) Ore ring (see also [271 Prop. 5.9]). 
Thus, we can form the field of (left or right) fractions K(x\ a, ■&). See for example [28] for an 
extensive introduction on how to properly define the various arithmetic operations for such 
fractions. 

If K is commutative, a = id and $ = are the identity and the zero map, respectively, 
then K[x; id, 0] is just the usual polynomial ring if [a;]. Two other important special cases of 
Ore polynomials are linear differential and delay operators, which we discuss in the following. 

Example 1 Let K be the field of meromorphic functions over the real line (see e.g. f29[ p. 

42])). 

1. Assume first that i? = and that a = 8 is the time delay operator which is defined 
by Sf(t) = f{t — t) for all f G K where t > is a fixed real number. In an abuse of 
notation, we will denote the ring K[x;8, 0] just by K[S\, i. e., we identify 8 with the 
Ore variable x. The commutation rule for K[8] is Sf(t) = f(t — t)6 for all fit) £ K. 
We call it the ring of time delay operators. 

2. Assume now, that a = id. Let t? = d be the usual derivation in the sense of calculus. 
Then, the ring K[x;\d,d\ — which we will just write as K[d] — has the commutation 
rule df(t) = f(t)d+ fit) for all fit) € K . This is the ring of differential operators. 

Since the maps d and 8 commute, we may extend d to K [5] by setting diS) = 0. Thus, 
the ring O = K [S, d] — or in more complete notation K [x; S, 0][y, id, d] — is well-defined and 
has the commutation rules 

Sf(t) = f(t-T)S, df(t) = f(t)d + f(t) and dS = Sd 

where fit) € K. We call O the ring of time-delay differential operators. The commutation 
of d and 8 means that this ring is an Ore algebra in the sense of [301 Def. 1.2]. Using the 
formulas in [28l Thm. 13], we may extend the action of d on the fractions in K (S). Thus, 
also the ring O = K (8) [d] is a well-defined Ore polynomial ring. 

Since 8 is an automorphism, by |25[ Sect. 7.3] we may form the ring of formal Laurent 
series K((8J) in 8 with coefficients in K. That is, K((8)) consists of elements of the form 
12j>N fj S3 where N e Z and fj e K for all j > N. Moreover, by [25j Prop. 7.3.7] we may 
embed K (8) into K (((5)). 
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2.2 Signal space 



The signal space has to be contained in the domain of the previously defined operators. 
Such a space is not unique in general and we may choose it in accordance to our control 
objectives. A signal space suitable for motion planning tasks should satisfy at least the 
following requirements. It should: 

1. be closed under the action of any polynomial of differentiation (d) and delay (5) 

2. be closed under advances defined by inverses of ^-polynomials, 

3. contain C°° functions with compact support in R. 
Therefore we introduce 

S = {/: R -> R | BEf C R discrete : / £ C°°(R \ E f ,R) 

and 3t £ K Vi £ K : * < t => /(*) = 0}. 

which is an adaptation of a space introduced in [3ll Sec. 3]. Indeed, S fullfills the first and 
third property and it is even an O-module. 

As far as property 2 is concerned: if it £ K[S] with ir ^ 0, let us compute ir^ 1 f for f £ S. 
To this end we can develop 7r _1 into a formal Laurent series n^ 1 — ^2j>N where iVeZ 
and aj £ K for all j > N (see |A| . Then for any / £ S and t e IR we have 

j>iV j>Af 

Since f(t — jr) = for jr > t — to all but finitely many of the summands on the right 
hand side vanish at any given t. We can check that 7r -1 / is again in S. Thus, property 2 
is satisfied. Moreover, one can show that S is an O-module. Note that changing O to O in 
conjunction with S will yield significant simplifications in the sequel. 

The example discussed in Section [5] also shows that S is relevant for motion planning 
tasks. Note in addition that S contains splines and step functions. 

3 Matrices &; hyper-regularity 

We model systems of linear time-delay differential equations using matrices of operators. 
The set of all n x m matrices with entries in O is denoted by O nxm . Square matrices 
in O nxn which possess a two-sided inverse that is also in O nxn are called unimodular. 
The set of all unimodular matrices of O nxn is denoted by G1„(C). We write 1„ for the 
ft x ft identity matrix and nxm for the n x m zero matrix. In both cases we will omit the 
indices when they are obvious from the context. We use lxm for the set of row vectors 
of length Til and O n for the set of column vectors of length n. Given a matrix M £ nxm 
we denote its O-row space by lxn M. For a matrix M — (M it j) £ O nxm we define 
degg M = max{deg a Mij | i — 1, . . . , n, j — 1, . . . , m}. This will also be applied to row or 
column vectors regarding them as 1 x m- or n x 1-matrices, respectively. We write Mj * for 
the i-th row of M and M* t j for the j-th column of M where i = 1, . . . , n and j = 1, . . . , m. 
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Let M G Qnxm^ g mce K(S) is a division ring, O is a principal ideal domain. This 
means we can apply |32l Thm. 8.1.1] in order to find unimodular matrices S € Gl n (C) and 
T G Gl m (0) such that SMT is in Smith- Jacob son form. This is a diagonal form with the 
additional property that each diagonal element is a total divisor of the next one. 

Definition 1 (Hyper-regularity |5j) A matrix M G O nxm is called hyper-regular if the 
diagonal elements of its Smith- Jacobson form are all 1, i. e., if there are unimodular matrices 
S G Gl n (0) and T G Gl m (0) such that 



'n—mxm 



ri > 771 and SMT = . 
or n<m and SMT = (1„, 0„ XTO _„). 

Usually, the given method for obtaining a Smith- Jacobson form involves computing a 
diagonal form by repeatedly taking greatest common (left or right) divisors as a first step. 
Since the performance of this approach is difficult to predict and most likely exponential 
due to degree growth, we will give an alternate, and computationally more efficient, charac- 
terisation of hyper-regularity below. 

Proposition 1 (|19j) (i) A matrix M G O nxm , with n < ra, is hyper-regular if, and 
only if, it possesses a right-inverse, i.e. 3 T G Gl m (0) such that MT = (l n , nx ( TO _ n )). 

(ii) A matrix M G Q nXm } with n > m, is hyper-regular if and only if it possesses a 

left-inverse, i.e. 3 S G G1„(0) such that SM = ( ^ m 

V"(n— m) x m 

Proof. We only prove (i). The proof of (ii) follows the same lines and is left to the reader. 
Let M G O nxm be given with n < m. M is hyper-regular if, and only if, there are matrices 
S G Gl n (0) and T G G\ m {0) such that SMT = (1„, 0„ xm _„). Thus, using the identity 



(i„, o) = s- 1 (1„, 0) 

we get 



S 
l m _ 



(1,„ 0)= S ->( S MT)g 1 „°_J= M ( r (o ," 
which proves that M is hyper-regular if, and only if, it has ^T (q ^ ^ ) ) ;is ri^lil-invcrsf. 



We will apply the technique of row-reduction for the computation of left- or right-inverses. 
A matrix M G Q nxm j s called row-reduced (or row-proper) if for all row vectors v € O 
the so-called predictable degree property 

deg a vM = max{deg a v + deg a M^* | j = 1, . . . ,n} 

holds. Note, that this definition differs from the usual one as given in 33, Sect. 2] or [34], 
Sect. 2.2], but is shown to be equivalent in [33l Lem. A.l (a)] where one can easily check that 
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the proof also works for division rings instead of fields. Similarly, the algorithm outlined 
in the proof of [33J Thm. 2.2] can be easily transferred to division rings, and we obtain 
that for every matrix M E Q nxm there exists a matrix S € G1„(0) such that the non-zero 
rows of SM form a row-reduced submatrix. Using the results in }35j row- reduction is of 
low polynomial complexity in the size of M and its degree. Directly from the definition we 
obtain that the rows of a row-reduced matrix must be linearly independent. 

Row-reducedness is connected to the Popov normal form (see, e.g., [351 Def. 2.2]) which 
is essentially a row-reduced matrix with additional properties to make it unique. One may 
prove that for each matrix in O nxm there exists exactly one matrix Popov form having the 
same row space. Also, row-reduction may be regarded as a special case of Grobner basis 
computation — see [55] , 

We consider now the case n < m. Let M e O nxm , and let S € G1„(0) be such that 



SM = 




where M S O is row-reduced and k is the (left) row-rank of M by [33J Thm. A. 2]. 
Assume first that M is hyper-regular. As discussed above, this means that all unit vectors 
of lxm are in the row-space of M. Thus, we must have k = m and by [531 Lem. A.l (c)] 
all rows of M have 9-degree or below. Thus, M £ K(5) m m . Since all rows of M are 
linearly independent, we conclude that M has maximal (left) row rank. 

Conversely, if row-reduction of M yields a matrix of degree and (left) row-rank m, then 
clearly M is hyper-regular. 

There is also the analogue concept of column-reduction. Each matrix may be brought into 
column-reduced form (up to zero rows) by right multiplication with a unimodular matrix. 
All the results cited above hold with the appropriate changes. In total, we have proved the 
following lemma. 

Lemma 1 ( [19j ) A matrix M 6 Q nxm j s hyper-regular if and only if 

1. n > m and row-reduction yields a matrix of d -degree and left row-rank m 

2. or, n < to and column-reduction yields a matrix of d -degree and right column-rank n. 



4 Systems 

Consider systems of the kind 

Ax = Bu, (1) 

with the pseudo state x of dimension n, input u of dimension m and the matrices A £ j^ nxn 
and B £ lZ nxm , where 72. is a ring. For the analysis of such systems, two important objects 
are considered (see e.g. [371 EH1 EE HOI III LI31 III HS1 LH] ) : 

• its behavior B = kei(A, —B), where the kernel is taken with respect to the chosen 
signal spac^] where the components of the variables x and u are supposed to live 

Hhis space is not uniquely denned and may be chosen according to the specific application one is interested 
in, such as motion planning, tracking, etc. See section|2| 
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• and its system module M = K lx ( n+m *> /K lxn (A, -B), where TZ lxp is the set of row 
vectors of length p, for every p £ N, with components in TZ and where TZ lxn (A, —B) 
is the module generated by the rows of the matrix (A, —B), 

Linear time-invariant differential systems, i. e. without delay, the ring TZ being chosen as 
R[d] and being commutative, are shown to be differentially flat if, and only if, their system 
module M. over TZ is free, and a flat output is, by definition, a basis of the free module 
M. (see e.g. [H[S]). This property may serve as a flatness definition for linear differential- 
delay systems ([7J [2D [B3])- Nevertheless, only few systems have a free system module. 
Moreover, since the considered base ring TZ = O = K[5, d] is not a principal ideal domain, 
freeness is in general different e.g. from torsion frccncss or projective freeness (see e.g. [7]), 
with important, but otherwise unclear, pratical consequences on the control possibilities 
of the system. However, for linear time-invariant differential-delay systems, namely with 
TZ = M[8,d], in virtue of the localization property of a commutative Ore algebra, torsion 
freeness of M. is equivalent to the existence of a so-called liberation polynomial tt =/= 0, 
7r £ R[<$], such that n^ 1 M. is free. This property is called 7r-freeness (see again [7]) and a 
basis of the free module 7r _1 M. is called a 7r-flat output. It can be interpreted as follows: 
if the system module, finitely generated by m input variables, is torsion free, it admits 
a basis of the form ir~ 1 {yi, . . . , y m } where yi £ TZ, i — l,...,m. Consequently, every 
system variable can be expressed as a combination of derivatives, delays and advances of 
(j/i, . . . , y m )- Note that the action of 7r _1 on a variable z may be interpreted as a formal 
power series in 5~ k z = z(- + kr) for k > —N, for some N £ N, and thus a combination of 
advances, or predictions, of z (see Section [2]). 

To extend this approach to time- varying differential-delay systems without restricting the 
time dependence of the system coefficients to be polynomial, as in [151 H2] where effective 
Grobner bases techniques are used, many new difficulties appear, and in particular the 
fact that the set made of the powers of a given polynomial is not in general an Ore set. 
Therefore, we are lead to propose a different definition of n-flatness, based on an extension 
of the "practical" flatness definition, in the spirit of [211 [S], saying that all the system 
variables are expressible in terms of a flat output, its derivatives, delays and advances in 
finite number. Moreover we introduce the definition of 7r-fc-flatness, which is an analogue, in 
our differential-delay context, to the notion of fc-flatness (see [2U BH 113 SI] ) ■ This viewpoint 
is developed in the next two subsections. 



4.1 Recalls on the framework 



From now on, we take System (JX) with TZ = O = K[d, 6] as in Section |2.1[ Recall that 5, 



the delay operator, is defined by Sf(t) = f(t — r) for all t £ R, where r is a given positive 
real number and f £ K, and that d = 4z is the ordinary time derivative operator. The 
ground field K is the field of meromorphic functions over the real line and the notation 



O = K(S)[d] is defined in Section 2.1 Furthermore, we will consider only the signal space 
S from Section HPI 

We assume that the matrix (A, —B) has full (left) row rank. 

In this case, System (JTJ is a differential-delay system, possibly with time-varying coeffi- 
cients, the coefficients being meromorphic functions of time. An example of such system is 
provided by the following: 
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Example 2 Consider the system 

x 1 {t) = a{t){x 2 {t-l)-x 2 {t-2)) 
x 2 {t) = u(t-l), 

where the pseudo-state x = (x±, x 2 ) belongs to S 2 , the control u belongs to S , and the function 
a is meromorphic with respect to time. Here, the matrices A and B are given by 

and, as announced, they are matrices over the ring O — K[S,d]. 
4.2 7r-flatness 

One usual way to consider System (JlJ is to bring all the variables in one side, i.e. to consider 
the system 

(A,-B) (Y)=0 



or, in other words, 



with F = (A, —B) and £ = 



Ft = (2) 



Definition 2 (7r-flatness) The system ^ is called 7r-flat if there exist tt G K[8] and ma- 
trices P e 0™ x («+ m ) and Q G {n+m ^ m such that 

iT^QS" 1 = B = kerF and 7T _1 P 7t _1 Q = l m . 

Equivalently, there exist a polynomial it € K[6] and matrices P and Q of suitable di- 
mensions over the ring O, such that £ = ir~ 1 Qy 1 with y — ir~ 1 P(; for all £ € S n+m , where 
y € <S m is a 7r-flat output. In other words, 7r-flatness means that all the system variables 
can be expressed as linear combinations of y and a finite number of its delays, advances 
and derivatives, in an invertible way, namely the matrix ir~ 1 P admits the matrix it~ 1 Q as 
right-pseudo-inverse. 

If 7r G K, then the system is simply called flat and y a flat output. 

This definition is thus an extension of the one proposed by [5J [3] for finite-dimensional 
nonlinear ordinary differential systems or by |21j for time- invariant delay-differential systems. 
Other definitions for time-varying linear differential-delay systems with polynomial time 
dependence are proposed by Chyzak, Quadrat and Robertz [To] . 

Remark 1 Note that Definition^is equivalent to the existence of matrices P G g} mx (" +7n ) 
andQ G O v ' such that QS m = B and PQ = l m . We recover tt by computation of 
the left common denominator of P and Q \21\ Prop. 5.3]. Therefore all computations can 
be done over O. 
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Wc have the following proposition: 

Proposition 2 Assume that y is a ix-flat output of system Q). Let us set y — Tz with 
T G Gl m (C). There exists a polynomial k G K[6] such that z is a K-flat output of system 

Proof. By Remark [TJ it is sufficient to prove that T~ 1 P and QT are related to the 
transformed flat output z. Obviously, we have T~ 1 PQT = l m and moreover QTS m = 
QS m = B. With k being a common denominator of T~ 1 P and QT, the claim follows. ■ 



Remark 2 To every -K-flat system there obviously corresponds a polynomial t:q G K[S] of 
minimal degree such that the system is TTQ-flat. 

To characterize 7r-flat systems, we introduce the following definition: 

Definition 3 We call M. = Q n+m / 'Q n p the extended system module. 

Theorem 1 We have the following equivalences: 

(i) System ^ is n-flat; 

(ii) The extended system module M. is free; 
(Hi) The matrix F is hyper-regular over O . 

Remark 3 Note that in the linear time-invariant case Definition [£| is equivalent to the 
definitions of P7| {TBjj. 

For the proof we need the following lemma: 

Lemma 2 Let A be an arbitrary ring and M € A pxq , and let S € Gl p (A) and T £ Gl g (A) 
be unimodular. Then A lx VA lxp M = A lx i/A 1x p(SMT) as left A-modules. 

Proof. The map ^ = !jh> vT is an automorphism of A lx? because T is unimodular. The 
kernel of the composition 7r o <p of the projection 7r: A lx<? — > A lxq /A lxp SMT with ip is 
ip- 1 (A lx PSMT) = A lxp SM = A lx PM where the second identity follows from the unimod- 
ularity of S. Since ir o ip is also surjective, by the first isomorphism theorem for modules 
(see, e.g., [13 Theorem 1.17]) there is an isomorphism A lx «/A lxp M ->• A lx * / A lxp (S MT) . 



Proof, [of Theorem [I] We prove first that ^ is equivalent to ( [iiij ) . If the system F£ = 

n - -rmx(n+m) - -^(n+m)xm . 

is 7r-fiat, there are matrices F G C and Q € C with common denominator 

7r G K[5] \ {0} such that QS m = B = ker F and P Q = l m . That means that Q is surjective 
and that the following diagramme is exact. 

S n+m > imF > 

^ ^ 

E 
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Since PQ = 1, by the splitting lemma (see, e. g., [J51 Prop. 3.2]) one can show (see [13]) that 

there is a matrix E € O such that FE = l n where n is the (left) rank of imF because 
the rows of F are linearly independent by assumption. Therefore F is hyper- regular. 

Conversely, let F be hyper-regular. Using Lemma[l]we may thus compute W € Gl„ +m (0) 
such that FW = (1„, 0). If we let Q be the last m columns of W and P the last m rows of 
TU _1 , then we have FQ — and PQ — l m . Extracting a common denominator it of P and 
Q, we have proved that Ft; — is 7r-flat. 

Next we show that ( |iii| ) is equivalent to (JTTJ) . Assume first that F is hyper-regular. Since 
the unit vectors are in the column space of F, by column-reduction we obtain an invertible 
matrix T € Gl n + m (0) such that FT = (1„, 0). By Lemmaj^J this means that 
-~lx(n+ro) -~lx(n+m) -^lx(n+m) 

-T-T _ U ^ U ^ _U ^ -pylXTO 

1X "F " lxn FT ~ O 1X "(l„,0) " 

which is free. 

Conversely, assume that M. is free. Using, e.g., the method in 32, Thm. 8.1.1] we may 
obtain unimodular matrices S £ Gl n (0) and T € Gl„+ m (0) such that SFT = (A, 0) where 
A = diag(a 1; . . . , a n ) G O is a diagonal matrix. By Lemma^l 

-=rlx(n+m) -=rlx(ra+m) -^rlx(n+ro) n ^ 

1X "F 1X "^T O lx "(A,0) 52 Oo," 

Since this module is free by assumption we conclude that all a,j must be units, i. e., we may 
assume w. 1. o.g. that A = 1„. By Definition [l] F is hyper-regular. ■ 

Algorithm 1 (Computation of a 7r-flat output) 
Input: A matrix F <E Q nx (. n + m ) representing the system 

Output: An Ore polynomial tt £ K[8] together with matrices P € Qmx(n+m) an( ^ q g 
as in Definition^ or FAIL if such matrices do not exist. 

Procedure: 

1. Use column-reduction to check whether F is hyper-regular. If not, then return 
FAIL. 

2. Else, let W € Gl n+m (0) be such that 

FW = (1„ 0„ xm ) . 

3. Let 

Q = t(H and P±(0 mxn l m )W-\ 

4- Let 7r 6 K[S] be a common denominator of P and Q. 

5. Set P = nP e mx («+"») and Q = ttQ € C>(™+ m ) x "\ 

6. Return tt, P and Q. 

Remark 4 In the above algorithm we exploit the fact that W~~ x can be computed at the 
same time as W by inverting the elementary actions that compose W . 
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4.3 7r-0-flatness 

In contrast to the considerations of the previous subsection, it is sometimes necessary to 
keep the state and input variables separate: in the theory of linear time-invariant systems, 
controllability is equivalent to the existence of Brunovsky's canonical form (see e.g. [UJ EE]); 
the interpretation of some of their states as flat output (see e.g. 0) shows that flat outputs 
do not need to depend on u, i.e. there exist P £ R mx " and Q £ R" xm sucn that y = Px, 
x = Qy and PQ = l m . This property is called 0- flatness. The fact that u can be expressed 
as a function of y is, in this case, an immediate consequence of the system equation with 
x = Qy. More generally: 



Definition 4 We say that a system is n-k-flat, with k > 1, if and only if there exists a 
TT-flat output y such that the maximal degree with respect to d of the matrix Pi 



>n—m 

1 i n 



equal to h — 1 . 

We set k = 0, by convention, if P |^"~ m j — o, i. e. y does not depend on u. 

Note that 7r-0-flatness is equivalent to the existence of (n, P, Q) as in Definition [5] such 
that P = (Pi m , ro ) with Pi £ O mxn and tt^Pitt^Qi = l m where Q x = (l„ 0„, m ) Q. 

A linear flat system is not necessarily 0-flat, as shown by the following elementary ex- 
ample: 

Example 3 Consider the system 

(1 -d) 



x 
u 

or 

X = u. 



It can be easily seen that y — u is a flat output. Hence the system is 1-flat. 

Remark 5 Note that elementary considerations show that 1-flatness is not preserved under 
the action of the group of unimodular matrices. Indeed, setting (^/j — ^ J * n the 

previous example, we get the system (—9 l) (^ij ~ or ~ u ' ' > which admits y' = x' 

as flat output and is thus obviously 0-flat. However, if we restrict the transformation group 
to be G\ n (0)®G\ m (0), a group preserving the control variables, then 1-flatness is preserved. 

Example [3] and Remark [5] thus show that obtaining a characterization of 7r-0-fiat systems 
is of interest. 

Lemma 3 (Elimination) If B in is hyper-regular, the system can be decomposed ac- 
cording to 

( 3^)x = MAx = MBu = ( n lm ) u. (3) 

F) \0( n _ m )x m / 

where F £ £>(»-™)x« and ip £ K[S] with (f^O. 
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Remark 6 Let F £ O . Let tp G K [6] be its common denominator. We have F — <p l F 
with F G O pxq . Thus, kerF = kerP. 

Theorem 2 If B is hyper-regular, we have the following equivalences: 
(i) The control system is TT-O-flat; 

(ii) The extended system module Q lxn JQ lxtyn p j s free, with F defined in d3l; 
fraj P is hyper-regular over O. 

Proof. Let F be hyper-regular over O. From representation ^ combined with Theorem [l] 
there exists Pi £ O™ "_and Q\ G O m such that FQi = and P\Qi = l m . Thus, 

P = (A, 0) and Q 4 ^ J satisfy PQ = = 1« and 

M (A -P) Q = M (AQ X - BRQi) = ( ^ F ) Qi - (Y) = °- ( 4 ) 

Equation @ implies that Q5 m C B. To prove the other inclusion, let j = ( ^ ) £ B. Then 
analogously 

= M(A -B)7=( (p -%)7r-{^)l, 

Since ker</? _1 P = kerP, we obtain 71 G kerP and P71 — 72. Consequently, there exists 
C G S m such that 71 = QiC an d thus also 72 = RQiC Hence, 7 G QS m and the other 
inclusion also holds. Thus, extracting a common denominator tt such that P = tt~ 1 P and 
Q = tt~ 1 Q where P G O mxn and Q G 0" xm , we see that tt, P and Q fulfill Definition g 
Thus, (iii) implies (i). 

Conversely, if (i) holds with tt, P and Q being as in Definition BJ then with Q\ — 
(1„,0)Q and Pi = P( 1 "), by the same calculation Q we obtain (p~Fn~ 1 Qi — and 
n~ 1 PiTr~ 1 Qi = l m . Thus, Fx = is tt- flat. Thus, by Theorem [l] we have (iii). ■ 

If the matrix B of system is hyper-regular, we can directly check for the existence of 
a 7r-0-flat output, as shown in Theorem [2j A suitable algorithm is as follows: 

Algorithm 2 (Computation of a 7r-0-flat output) 

Input: Matrices A G Q nxn an d q g Qnxm w ^ £j hyper-regular, representing 

Output: If system is TT-O-flat, the polynomial tt G K[S] and the triple (P, Q, R) of the 
matrices P G O mxn , Q G O nxm and R G O mxm from Definition^ 

Else, if the system defined by A and B is not TT-O-flat, then FAIL. 
Procedure: 

lm 

(n—m) xm 



1. Compute M G G\ n {0), s.t. MB = . 
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2. Write 

R 



MA = 



ip- L F 



where Re O and F £ Q{n-m)xn ^ 
3. If Algorithm^ applied to F returns tti, P\ and Q%, then 

(a) Set P = (ir- 1 P 1 ,0) and 

n - ( n ^ Ql ' 

(b) Let 7r be a common denominator of P and Q and set P — ttP 6 O r 
Q = irQe O nxm . 

(c) Return 7r , P and Q. 
4- Else, return FAIL 

If the matrix B of system ([!]) is not hyper-regular, we can check for 7r-fiatness by applying 
Algorithm [T] 

Remark 7 Note that Algorithms^and^do not necessarily yield a polynomial ir of minimal 
degree in S. Especially, it is not guaranteed, that we obtain ir £ K if the system is flat. 

5 Example 

To illustrate the results of the preceding sections and the usefulness of the concept for the 
feedforward controller design for linear time-delay differential systems, we demonstrate all 
steps for the system of Example [2j Note that all the necessary computations can be done 
with a package for the computer algebra system Maple, which has been presented in |19j . 
The package can be obtained from the first author upon request. 

The matrix B is hyper-regular and following Algorithm [2] we compute M such that 

MB = Q . We get 

s- 1 " 

1 



M = 
Thus: 



We now follow Algorithm [l] and compute W such that FW = (l 0) 

w = {(s 2 -V" -(<P -V 1 ^) ' 

The inverse is 

'd a(S 2 -S) 



to 



1 
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Subsuming the remaining computations of Algorithm [T] we directly obtain the following 
matrices in step 3 a) of Algorithm [2j namely 

P=((0 l)^" 1 0) = (1 0), (5) 

and 

-^-yi 9 V (6) 
v ' y j \-(?-*r l (&-$d)j 

We readily get tt = S 3 — 5 2 . 

From ([5]), a 7r-0-flat output is given by y = X\. This can also be deduced from the first en- 
try in Q. From (JgJ we deduce that x 2 — -(5 2 — S)^ 1 and u = — (5 3 — S 2 ) 1 (Ay — ^y). 

We now use these formulas to generate feedforward trajectories corresponding to the 
following motion planning problem: we are looking for a trajectory of {xi,X2,u) starting 
from x(0) = (xi(0), £2(0)) = (0 0) and arriving at the final state x(2) = (1 0) at time t = 2, 
while X\ has to be equal to zero before t = and to 1 after t = 2. 

In order to compute explicit feedforward state and input trajectories without integrating 
the system differential-delay equations, we compute a reference trajectory for the flat output 
y, which is not constrained by any differential-delay equation thanks to the freeness of the 

extended system module, and then deduce the state and input trajectories by = Qy 

with Q defined by (JsJ) . We take a(t) = t + 3 and thus a(t) = 1. Furthermore, we take t = 1, 
i.e. Sa(t) = a(t — 1). 

In order to obtain an explicit expression for x 2 — ~{S 2 — ( v)> we compute the 
series expansion 

(s 2 -sy 1 = J2 -v- 
j>-i 

Taking into account that x\ is required to be constant for t < and t > 2, the desired 
trajectory yd for y = X\ will be constant at these points of time and thus yd = jjd = for 
t < and t > 2. This yields 

X 2 ,d = ( P J = X2 > d (*)X[-1,0I (t) + (*2,d(*) + ^2,d(< + l))X[0,oo[(* + 1), (7) 



where X 2 , d (t) = a(t \ t]) y<i{t -\t\). 
We make the ansatz 



y d (t) = (a + ait + a 2 t 2 + a 3 t 3 + a 4 t 4 + a 5 t 5 )x[ ,2[(t) + Xp,oo[(*) 

the coefficients dj, i = 0, . . . , 5 being such that y(0) = xi(0) = 0, y{2) = x\{2) = 1, y(0) = 0, 
2/(2) = 0, y(0) = and y(2) = 0. We readily get: 

/ 45 35 3 \ 

= + ^ - X[0,2[(*) + X[2,oo[(*)- (8) 
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We now insert this expression in ([7]) to obtain the corresponding feedforward control. Since 

7T- 1 = (S'-S 2 )^ = E &3 ( 9 ) 

J>-2 

we get 

u d {t) = U d {t) X[ -2,-i[(t) + (U d (t) + U d (t + l))x[-i,oo[(0. ( 10 ) 

where U d {t) = - a(t ^ tm ij d (t - [fj) + a 2 (t-Lt]) &*(* ~ L*J ) •> which completes the solution of 
this motion planning problem. The resulting trajectories are shown in Figure [1] Though the 
trajectory for y is relatively smooth, it is generated by an input u which is quite irregular and 
complicated. We also notice that u has to start 2t before the beginning of the y-trajectory 
in order to satisfy the above requirements, which corresponds to the order of the Laurent 
series associated to 7r _1 . 






Figure 1: Graphs for X\ = y, x-i and u 



6 Conclusion 

In this paper we have introduced a new definition of 7r-flatness and the important particular 
case of 7r-0- flatness. These notions have been characterized by means of the properties of ma- 
trices over Ore polynomial rings, such as hyper-regularity or row- and column-reducedness. 
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This characterization has yielded an efficient and simple algorithm to test the existence of 
7r- and 7T-0-fl.at outputs and compute them, via an algorithm that is implemented in the 
computer algebra system Maple |19j . We have illustrated the pertinence of this approach 
by an academic example. 

Note that the description of all possible 7r-0-flat outputs can be achieved; it will be the 
subject of a forthcoming paper. This will lead to a method to determine the minimal order 
of the polynomial n. 

Another potential direction to extend this approach may be to investigate other possible 
signal spaces. 

Finally, the question of increasing the efficiency of the proposed algorithm by using 
quotient free computations may also be addressed. 
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A Inversion of a ^-Polynomial 

Let 7r = a n S n + . . . + a k 5 k e K[S] where k < n and a k , . . . , a n G K. Assume that a k ^ 0, 
i.e. that the order of 7r is k. The inverse ir^ 1 in K((5)) may be computed in the following 
way: First, using the fact that 

w- 1 =5- k (a n 5 n - k + ... + a k )- 1 

we may assume w. 1. o. g. that k = 0. We make an ansatz J2j>o CjS^ for n^ 1 and compute 

n min{^,n} 

l = 7r^ C ^=^E a ^fe)^=E( E ^(ce-i))^. 
j>0 i=0 j>0 t>o i=0 

Comparing coefficients and using oo ^ 0, we can compute 

1 = a c c = a^ 1 

and 

min{^,n} min{^,n} — 1 

= ^ a i 5 t {c t -i) ce = -% 1 ^ aiS % (ce-i) 

i=0 i=0 

for all I > 1. Note, that the left hand side depends only on those Cj which are already 
computed. 
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